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ABSTRACT 

A  DIRECT  METHOD  FOR  THE  MEASUREMENT 
OF  CONTROL  ROD  WORTH   IN   A  NUCLEAR  REACTOR 

By 

William  Arthur  Emslie 


It  is  often  necessary  to  perform  control  rod  calibration 
in  a  nuclear  reactor.   Such  calibration  requires  the  measurement 
of  the  reactivity  produced  when  discrete  sections  of  a  control  rod 
are  withdrawn  from  a  critical  reactor  pile.   From  these  data  a 
total  reactivity  worth  curve  can  be  complied  for  the  length  of  the 
rod.   The  subject  of  this  thesis  is  the  development  of  a  means  by 
which  reactivity  worth  can  be  measured. 

A  solution  to  this  measurement  problem  has  been  developed 

using  the  current  output  of  a  compensated  ionization  chamber 

(adjacent  to  the  reactor)  as  a  signal  source.   When  a  section  of 

control  rod  is  withdrawn  from  a  critical  reactor,  the  current  from 

the  ion  chamber  is 

t/8 
ls(t)  =  I0e   P 

where      t      is    time,    9        is    the   reactor   stable    period   and      I_       is 

p  0 

the    initial    ion   chamber   current. 

The   stable   period    is   proportional    to    the    amount   of    reactivity 
produced   since 
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9p        k       -  p 
ex 

9^   is  the  mean  neutron  lifetime,  k     is  the  excess  multiplication 
0  ex 

factor  and   p   is  the  reactivity.   This  approximate  relationship 
between  9   and   p  makes  it  possible  to  obtain  the  reactivity 
from  the  output  current  of  the  ion  chamber  since 

t/e  pt/en 

Actual  measurement  of  p  is  accomplished  with  a  circuit 
composed  of  a  logarithmic  amplifier  followed  by  a  differentiator 
and  filter.   The  transfer  function  of  this  circuit  is 


(-R^s) 
V  As)    -  V,(s) 


outv  '     lv  '  (R1C1s+l)(R2C2s+l)   (RfC  s+1) 

where     V.(s)       is    the    frequency   domain   expression   for     V    (t)       (the 
output  voltage    of   the    logarithmic   amplifier)    and    is    equal    to 

v    ,   ,         -A    log10(IQ/lk)         A  log1Qe 

V  (a)  = —2- 

9  s 
P 

In  the  time  domain,  the  logarithm  of  the  ion  chamber  current 
produces  a  linear  ramp  voltage  v.(t)   where 


t/e 

i  (t)  I  e   p  I 

v1(t)  =  -A  log1()  -A—  =  -A  log   — —  =  -A  log1()  -  -  ^  log1()e 

k  k  k    p 

This  ramp  voltage  is  then  differentiated  and  filtered  to 
obtain  a  d-c  voltage  proportional  to  the  inverse  of  p.  The  re- 
sulting voltage   (V«)   is 
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v2  "  "K  37  (vi(t»  "  "K  &  ('A  logio  17  -  f  logioe> 

k         p 

KA   .  oKA   , 

=  r  logioe  a  \  logioe 

p  o 

V„  may  be  calibrated  on  a  d-c  voltmeter  to  indicate  the 
corresponding  reactivity. 

A  circuit  to  synthesize  the  transfer  function  was  designed 
and  constructed.   Results  of  tests  made  on  circuit  performance  show 
the  total  accuracy  in  the  measurement  of   p   to  be  within  +0.5 
percent  of  the  meter  scale. 
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CHAPTER  I 
INTRODUCTION 

The  necessity  to  perform  calibration  of  the  control  rods 
in  a  nuclear  reactor  requires  a  means  by  which  to  measure  the 
reactivity  produced  by  withdrawal  of  discrete  sections  of  a  con- 
trol rod  from  a  critical  reactor  pile. 

Use  of  an  ionization  chamber  adjacent  to  the  reactor  pro- 
vides an  exponential  current   (i  (t))   from  which  the  reactivity 
may  be  derived.   When  a  section  of  control  rod  is  withdrawn  from 

a  critical  reactor   i  (t)   becomes 

s  ' 

t/e 

is(t)=IQe   P  (1.1) 

where   I_   is  the  initial  ion  chamber  current,  t   is  time  and   Q 
0  p 

is  the  reactor  stable  period. 

The  reactivity  may  be  obtained  from  equation  (1.1)  since 

e„^—  (i.2) 

p     p 

2 
where   9n   Is  the  mean  neutron  lifetime. 


The  reactor  stable  period  used  in  this  expression  is  the  time 
required  for  the  ion  chamber  current   (i  (t))   to  increase  by  a 
factor  of  e. 

2 

This  approximation  is  explained  and  justified  in  Section  2.1. 


The  relationship  between  p   and   i  (t)   expressed  in 

s 

Equations  (1.1)  and  (1.2)  makes  it  possible  to  obtain  p   from 
i  (t)  .   This  may  be  accomplished  by  means  of  a  circuit  which  first 
takes  the  logarithm  of   i  (t)   and  then  differentiates  and  filters 
the  result  to  yield  a  d-c  voltage  proportional  to   p. 

An  example  of  a  control  rod  calibration  curve  is  illustrated 
in  Figure  1.1.   The  curve  is  compiled  by  withdrawing  the  rod  a 
short  distance  and  measuring  the  corresponding  reactivity  produced. 

This  procedure  is  repeated  until  the  entire  length  of  the  rod  has 

3 
been  withdrawn. 

Chapter  II  of  this  study  details  the  theory  underlying 
reactivity  and  the  ionization  chamber.   The  effects  of  prompt  and 
delayed  neutrons  are  also  discussed. 

In  Chapter  III  the  reactivity  measurement  problem  is 
formulated  with  a  set  of  possible  solutions.   A  particular  solution 
method  from  these  alternatives  is  then  presented.   Analysis  of  the 
design,  mathematical  properties,   and  calibration  of  the  circuit  follows 
Further  design  analysis  is  presented  in  the  Appendices. 


3 

The  units  of  reactivity  as  expressed  in  Figure  1.1  are  in  dollars 

and  cents  where  1  dollar  =  100  cents  and  is  defined  as 

k 
1  dollar  =  -££  s-  £  (1.3) 

0.0. 

i   *i 

where   k    is  the  reactor  excess  multiplication  factor  and   ft.  is 

ex  v  >i 

a  normalizing   constant   equal   to   the    total    fraction   of  delayed  neutrons 
and  varies   with  different    reactors.      For   the  Michigan  State   University 
TRIGA  Mark   1  reactor     (3^=  0.0073      and   therefore   one   dollar's  worth 
of  reactivity    is    that   amount   required   to  produce   an  excess 
multiplication   factor  of     0.0073. 
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Figure   1.1      Calibration    curve    for    the    shim   rod    in    the  Michigan 
State   University  TRIGA  Mark  I    nuclear    reactor 


Chapter  IV  describes  a  prototype  instrument  which  was  built 
from  the  design  in  the  previous  chapter.   Test  results  and  circuit 
specifications  are  also  given. 

Conclusions  are  made  in  Chapter  V. 


CHAPTER  II 
REACTOR  TRANSIENT  BEHAVIOR 

4 
2.1   Reactivity,  Reactor  Period  and  Theoretical  Equations 

The  excess  reactivity  (hereafter  referred  to  as  reactivity) 
in  a  nuclear  reactor  is  defined  as  the  ratio  of  excess  to  effective 
multiplication  factors  and  is  given  by 

k 
eff 
In  the  above  formula  the  effective  multiplication  factor, 
k  f   is  equal  to  unity  for  a  critical  reactor.   The  excess  multi- 
plication factor  is 

k   =  k  _  -  1  (2.2) 

ex    eff 

and    is    equal    to   zero   for   a   critical   reactor. 

Substitution   of  Equation    (2.2)    into  Equation    (2.1)    yields 

keff    '    l 

P  -  -nr (2-3) 

eff 

Since    the   values   of      k    _,      used    in   this    analysis    are    to  be   no 

eff 


4 

The    content   of   this    section,    except    for    that  material   noted 

otherwise,    is    taken   from  M.M.   El-Wakil,   Nuclear  Power  Engineering 

(New  York:     McGraw-Hill  Book  Co.,   1962),   pp.   128-130. 


greater  than  1.03  ,  Equation  (2.3)  may  be  approximated  as 


Using  the  above  approximation  for  reactivity  it  is  possible 
to  relate   p   to  the  instantaneous  neutron  flux  within  the  reactor. 

The  equation  describing  the  flux  in  a  supercritical  reactor 
is 

t/e 

$  =  $Qe   P  (2.5) 

The  reactor  stable  period   (Q  )   may  be  written  as 

eo 
e    =  r~  (2.6) 

v  ex 

and   substitution  of  Equation    (2.4)    for      k  relates      p      to      Q 

ex  r  p 

where 

e    2=-"-  (1.2) 

p      p 

This  approximate  expression  for  9  may  be  substituted 
into  Equation  (2.5)  to  form  a  relationship  between  reactor  flux 
and  reactivity  such  that 

pt/e 

$  2=  $Qe    u  (2.7) 


Reference    to   the    rod   calibration   curve    in  Figure    1.1   shows    the 
maximum  reactivity   to   be   approximately   four  dollars.      This   value 
corresponds    to  an  effective  multiplication   factor    of   1.0292   as 
defined    in   the  Nuclear  Reactor  Operations    and  Training  Manuel, 
Reactor  Theory    (Michigan  State   University,    1969),    pp.    4-5. 


2.2  Relationship  Between  Reactivity  and  the  Reactor  Stable  Period 

As  Equation  (1.2)  illustrates,  the  approximate  relation- 
ship between  the  reactor  stable  period  and  the  reactivity  is  weighted 
by  the  mean  neutron  lifetime.   This  equation  is  characterized  by  a 
nonlinear  property  of   9   which  is  not  constant  but  rather  varies 
slightly  with   p.   Examination  of  the  exact  expression  for  reactivity 
reveals  the  source  of  this  nonlinear ity.   p  may  be  written  exactly  as 

90      V     Pi 
p  =  — r1 —  +  /  * (2.8) 

p  eff    i       l  p 
where   @.   is  the  fraction  of  delayed  fission  neutrons  in  the  ith 
group  and   \.   is  the  respective  decay  constant. 

For  large   Q  ,  Equation  (2.8)  may  be  approximated  as 

p   l    l      p 

This  approximation,  however,  is  only  valid  for   9  >  100   and  as 

P 

9   decreases  in  value  both  the  first  term  on  the  right  in  Equa- 
tion (2.8)  and  the  1  in  the  denominator  of  the  summation  no  longer 
remain  negligible.   The  net  effect  is  a  nonlinear  relationship 
between   9   and   p   as  illustrated  in  Figure  2.1. 

2.3  Effects  of  Prompt  and  Delayed  Neutrons  on  Reactor  Transient 
Behavior 

The  fission  process  in  the  Michigan  State  University  TRIGA 
Mark  I  nuclear  reactor  is  described  by  the  following  nuclear 


Samuel  Glasstone,  Nuclear  Reactor  Engineering  (New  York: 
D.  Van  Nostrand  Co.  Inc.,  1963),  pp.  244-245. 

7   Ibid.,  pp.  244-245. 
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Figure  2.1   Relationship  between  reactivity  and  reactor  stable 

period  for  the  Michigan  State  University  TRIGA  Mark  I 
nuclear  reactor  (Courtesy  of  Gulf  General  Atomic) 


equations 


92U235  +  Qnl  -   92U236  +  energy  (2.10) 


or 


235     1  1 

92U    +  0°  "*  2  fipS  +  enerSy  +  V  +  c0n  (2.11) 

Equation  (2.10)  occurs  about  twenty  percent  of  the  time.   The 

fission  products  (fips)  in  Equation  (2.11)  have  a  variety  of  atomic 

235 
weights  but  will  generally  weigh  about  half  as  much  as  the  U 

Heat  is  generated  by  the  energy  given  off  in  both  reactions  and  the 

energy  is  about  200  Mev  per  fission.  On  the  average  between  two  and 

three  (2  <,   c  ^3)  neutrons  are  given  off  per  reaction  in  Equation 

(2.11).   This  equation  is  critical  to  reactor  operation  since  the 

o 

neturons  produced  must  sustain  the  fission  chain  reaction. 

Neutrons  emitted  in  Equation  (2.11)  are  of  two  types: 
Prompt  neutrons  and  delayed  neutrons.   Prompt  neutrons  make  up 
about  99.27  percent  of  all  neutrons  emitted  in  the  TRIGA  reactor. 
The  total  mean  lifetime  of  these  prompt  neutrons  is  less  than  one 
millisecond. 

Delayed  neutrons  are  the  result  of  the  decay  of  certain 
fission  products  in  Equation  (2.11).   An  example  of  a  delayed 
neutron  emission  is  given  by  the  nuclear  equation 

53i137  _l  5/37   _^    5y36  (2-i2) 

->-*     oo — *       54  e    .   -  -  .  >      54  e 
22  sec        instantaneous 


8 

Training  Manuel,  Reactor  Theory,  p.  1. 

9  Ibid.,  p.  4. 


10 


137 
where   --I      is  one  °f  tne  fission  products.   The  22  second  delay 

in  the  neutron  emission  gives  it  a  half  life  of  appearance  of  22 

A        10 

seconds . 

Even  though  the  fraction  of  delayed  neutrons  is  small,  their 
long  apparent  lifetimes  make  the  average  lifetime  of  all  neutrons 
emitted  much  longer  than  the  lifetime  of  prompt  neutrons  alone. 
This  effect  is  vital  since  a  reactor  period  caused  by  prompt  neutrons 
alone  would  make  control  of  the  fission  process  impossible. 

The  effect  of  prompt  and  delayed  neutrons  on  the  change  in 
neutron  flux  can  be  described  as  follows:   In  a  critical  reactor 
withdrawal  of  a  control  rod  some  discrete  distance  will  produce  a 
transient  period  due  to  prompt  neutrons  and  a  stable  period  due  to 
the  combined  effect  of  prompt  and  delayed  neutrons.   The  effect  of 
the  transient  period  will  result  in  a  "bulge"  in  the  neutron  flux. 


This  effect  is  shown  graphically  in  Figure  2.2. 
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Figure  2.2  Combined  effect  of  prompt  and  delayed  neutrons  in 
a  supercritical  reactor 
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El-Wakil,  Nuclear  Power  Engineering,  p.  131 
Ibid.,  pp.  133-134. 
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Since   p   is  proportional  to  the  slope  of  the  curves  in 
Figure  2.1  it  is  necessary  to  let  the  transient  period  die  out 
prior  to  measurement  of   p. 

2.4  The  Compensated  Ion  Chamber 

It  has  been  stated  in  Equation  (2.7)  that  the  reactivity 
produced  when  a  discrete  portion  of  control  rod  is  withdrawn  from 
a  critical  reactor  pile  is  exponentially  proportional  to  the  change 
in  neutron  flux.   In  order  to  measure   p   it  is  therefore  necessary 
to  monitor  this  flux.   This  is  accomplished  in  the  TRIGA  Mark  I 
reactor  by  means  of  the  output  current  of  a  compensated  ion  chamber 
located  near  the  reactor  core. 

The  output  of  the  ion  chamber  is  a  very  low  current  which 
is  directly  proportional  to  the  neutron  flux  in  Equation  (2.5) 
such  that 

t/e 

l.(t)  =  IQe   P  (1.1) 

where   9   may  be  substituted  by  the  approximate  expression 

en^—  d.2) 

p   p 

as   previously  stated. 

Figure    2.3    is   a  simplified   diagram  of    the    compensated    ion 
chamber.      The   device    is   actually   composed   of  two   chambers   where 
R        is    the   equivalent    output    resistance.      One    chamber    is    lined 
with  boron  enriched  with  B        and    is   biased  at  +  580  volts.      From 
this   sub-chamber   a  current   proportional  to  neutron  flux   and  gamma 
flux    is    obtained.      The   other   sub-chamber   has    a  variable   bias    from 


12 


0  to  -37  volts  and  is  sensitive  to  gamma  flux  only.   If  the  gamma 
sensitive  current  is  subtracted  from  the  gamma,  neutron  current 

of  the  boron  coated  sub-chamber,  a  current  proportional  to  neutron 

12 
flux  alone  is  obtained. 
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Figure  2.3   Simplified  diagram  of  a  compensated  ion  chamber 

2    2 
At  low  neutron  flux  levels  ($  <  2.5  X  10  n/cm  /sec),  the 

current  contribution  from  the  gamma  sensitive  chamber  is  nearly 

equal  to  that  of  the  gamma,  neutron  sensitive  chamber.   It  is 

therefore  important  to  set  the  compensating  chamber  voltage  (V  ) 

to  the  proper  value  in  order  to  achieve  a  linear  response  due  to 

neutron  flux  alone.   This  leaves  room  for  error  in  the  linearity 

of   i  (t)   if  V    is  not  properly  adjusted.   Fortunately,  as  the 

neutron  flux  level  increases  the  current  contribution  to   i  (t) 

s 


12 


Training  Manuel,  Instrumentation,  p.  5. 
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due  to  gamma  flux  alone  becomes  negligible  and   i  (t)   approaches 
a  linear  dependance  on  neutron  flux  alone. 

Experimental  results  reveal  that  improper  adjustment  of 
the  compensating  voltage  at  low  reactor  flux  levels  (corresponding 
to  reactor  output  power  levels  below  ten  watts)  causes  the  output 
current  of  the  ion  chamber   (i  (t))   to  become  nonlinear.   Figure 
2.4  shows  the  maximum  nonlinear  effects  at  low  power  levels  when 
the  compensating  voltage  is  adjusted  to  zero  volts. 

Since  a  linear  relationship  between  ion  chamber  current 
(i  (t))   and  neutron  flux  (also  reactor  output  power  level)  is 
necessary,  it  is  best  to  rely  only  upon  output  currents  which 
correspond  to  reactor  power  levels  greater  than  ten  watts. 

Specifications  of  the  ion  chamber  are  given  in  Appendix  E. 

13 
Since  the  output  resistance  is  greater  than  10   ohms,  the  chamber 

may  be  considered  as  a  current  source. 
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Figure  2.4  Graph  of  ion  chamber  output  current  as  a  function  of 
output  power  level  for  the  Michigan  State  University 
TRIGA  Mark  I  nuclear  reactor 


CHAPTER  III 
DESIGN  OF  THE  REACTIVITY  MEASURING  DEVICE 

3.1  Formulation  of  the  Problem 

Having  presented  the  equations  defining  reactivity  and  its 
relationship  to  the  output  current  of  the  ion  chamber  it  now  re- 
mains to  develop  a  means  by  which  p   can  be  measured. 

Equation  (1.1)  expresses  the  relationship  between   9 

and  I   where 
s 

t/e 
is(t)  =  i0.     P  (i.D 

The  objective  is  to  process  this  current  so  as  to  produce  a  voltage 
proportional  to  the  reactor  stable  period.   If  this  can  be  accomplished 
effectively  then   9   may  be  related  to   p   as  in  Equation  (1.2) 
since 

9p^  (1'2) 

Because  this  relationship  is  slightly  nonlinear  as  the 
graph  in  Figure  2.2  illustrates,  the  value  of  p   corresponding 
to  a  given  9   may  be  taken  from  this  graph  and  used  to  calibrate 
the  reactivity  meter. 
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3.2     Possible  Methods    of   Solution 


13 


Given  the  desired  output   p   and  the  input   i  (t)   several 

methods  for  determining   p  will  now  be  considered. 

One  method  makes  use  of  the  circuit  in  Figure  3.1  where  a 

d-c  voltage  proportional  to   1/9    is  obtained  from   i  (t)   using 

p  s 

a  logarithmic  amplifier  followed  by  a  differentiator  and  filter. 


i    (t) 
s 

Logar i 

thmic 

vx(t) 

Diff 

sren- 

V2 
n 

Low-P 
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V 
01 

Amplifier 

tiator 

Filter 

fJ 

_L 

L 

Meter 


Figure  3.1   Block  diagram  of  a  circuit  which  provides  a  d-c 
voltage  proportional  to  the  inverse  of  the 
reactor  stable  period 

In  a  circuit  of  this  type  the  logarithm  of  the  ion  chamber 
current  is  first  taken  such  that 

vc>  -  -A  ^oV  =  -A  logioV  p  -  -A  iogiovt;  logioe  OA) 

The  voltage      v,  (t)       is    a    linear   ramp   function  whose    slope 

is    inversely   proportional    to   the   negative    inverse    of      Q    .      Dif- 

P 

ferentiation  of    this    ramp  voltage   yields    a   d-c  voltage     V        where 

~Kd  /    n  ~Kd    /    a    1  T      At    i  n  KA    1  /Q    OX 

v2  =  "dF  vi(t)  =  "It"  ("A  ^io^T  logioe)  =  T  logioe        (3'2) 

p  p 
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The  first  two  methods  presented  were  suggested  by  Dr.  L.  Giacoletto 

of  Michigan  State  University.   The  third  method  was  developed  by 
the  author. 
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The    result    is    a   d-c  voltage    inversely   proportional   to      8    .      This 

voltage    is    then  passed   through   a    low-pass    filter    to  eliminate    any 

low-frequency   noise    components   passed  by   the   differentiator. 

Finally,    due    to    the    slight   nonlinear   relationship   between 

ft        and      p      the  voltage     V  (Figure   3.1)    is    calibrated   on   a 

p  K  6  out      v      6  ' 

meter  via  Figure   3.6  to   indicate   exact   values   of   reactivity. 

A   second  method    for   obtaining   a   d-c  voltage   proportional 

to      9        involves    the   use    of   an    integrator   and   a  quarter-square 
P 

multiplier.      The    circuit    is   pictured    in  Figure   3.2. 


I 
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1 
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Figure   3.2     Circuit    for   obtaining   a   d-c  voltage   proportional 

to      9 
P 

In  this  circuit  the  input  current   (i  (t))   is  first  con- 

s 

verted  to  a  voltage 


t/e 

vl(t)  =  is(t)Rv  =  I0Rve    P 


(3.3) 


using   a  buffer   amplifier  whose    input   resistance    is   much   greater 


than     R    . 
v 


This  voltage  is  then  integrated  to  obtain  v„(t)   where 


i„r  9       t/e 


V^eIV^— n;P^    p-» 


(3.4) 
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The   quarter-square  multiplier    is   next   used   to  divide    out 
the   exponential    term    in     v    (t)      such    that 

-e 


v   (t)        -LRJ         t/e  , 

V  =  -1 =        °  v   P    (e        p    -    n    - 

out        v.(t)  RC         K  }  t/e 

T    R   e        P 
0   v 


RC 


(3.5) 


This  resultant  d-c  voltage  is  proportional  to   Q   and 

P 

can  therefore  be  calibrated  on  a  meter  using  the  graph  of  Figure 

3.6  to  represent  the  corresponding  values  of  reactivity. 

A  third  method  for  finding   9    involves  the  integration 

P 

of  a  voltage  proportional  to  the  change  of  the  input  current. 

Figure  3.3  is  a  simplified  diagram  of  the  circuit. 

v  (t) 
c 


Window 
Comparator 
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GO 


Reed 


Switch  |~Ly\M-- 


Figure  3.3   Circuit  for  obtaining  a  discrete  voltage  propor- 
tional to  the  slope  of   i  (t) 

s 

The    input    current    is    first   converted   to   a  voltage      v    (t) 

such    that 

t/e 


v    (t)    =    i    (t)R      =  I  R   e 
i  s  v  0  v 


(3.3) 


This  voltage  is  then  fed  into  a  window  comparator  whose  "GO"  out- 
put (Figure  3.3)  acts  as  a  current  sink  for  V  ^  v. (t)  £  V  and 
as  a  large  resistance  whenever  v,  (t)  <  V    or  v  (t)  >  V  .   The 

"GO"  condition  will  allow  the  current   i  (t)   to  flow  in  the  coil 

c 
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of  the  reed  switch  and  close  its  contacts.   When  this  occurs  the 

d-c  voltage  V    is  switched  to  the  input  of  the  integrator  and 

is  integrated  for  the  period  of  time  that   v..  (t)   remains  within 

the  bounds  of  V    and  V  .   This  period  of  time  is  proportional 
u        I 

to   9   and  is  also  proportional  to  the  final  integrator  output 
P 

voltage.   The  relationship  between  V     and   9   can  be  found 

out         p 

by  starting  with  the  output  voltage  integral 


-1  t 
V  t   =  5F  F  Vidt  (3-6) 

out    RC  d    I 


Th 


e  total  integration  time   (9  )   is  that  period  of  time  that 

V  <.   v.  (t)  £  V   and  may  be  calculated  by  substituting  V    and 

V  into  the  equation  for  v..  (t)   where 

W  -»*-Vo«  (3-7) 


and 


V8u>  "  Vu  ■  RvV  U   P  (3-8) 

9   and   9   are  the  times  for  which  v  (t)   is  equal  to  V   and 

V        respectively. 

u  r  J 

Equations    (3.7)    and    (3.8)    are   then  solved   for   the   time 

increment      9       where      9=9      -   9        and 
t  t  u         {, 

\  -  \  -  \  •  v*n  h;  "*n  v?  (3-9) 

Using  Equation  (3.9)  for  the  value  of   9    in  the  integral 
of  Equation  (3.6),  the  voltage  V     may  be  written  as 

vout  =  £  Vtn  fr  •*»  It?  <3'10) 

v  v  0      v  0 
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which  is  a  d-c  voltage  proportional  to  the  stable  reactor  period. 

As  in  the  first  two  methods  this  voltage  may  be  calibrated 
on  a  meter  using  the  nonlinear  graphical  relationship  of  Figure  3.6 
to  indicate  the  corresponding  value  of  reactivity. 

The  above  method  for  measurement  of  reactivity  differs 
from  the  preceeding  two  in  that  the  output  voltage  is  not  a  con- 
tinuous representation  of   9   but  rather  a  discrete  value  obtained 
after  integrating  a  d-c  voltage  over  a  discrete  time  period   9  • 

After  consideration  of  the  three  alternatives  the  method 
chosen  to  measure   p  was  the  first  one  described  using  a 
logarithmic  amplifier  followed  by  a  differentiator  and  filter 
(Figure  3.1). 

3.3  The  Circuit  and  Theory  of  Operation 

With  the  reactivity  measurement  solution  method  chosen, 

the  circuit  of  Figure  3.4  was  designed  to  produce  a  d-c  voltage 

proportional  to   9  . 

P 

The  transfer  function  of  the  circuit  is 

("R2C1S)  1 

Vout  =:Vl(s)  (R]C  s+1)  (R2C2s+l)  (RfCfs+l)  (3-11) 

where  V..  (s)   is  the  frequency  domain  expression  for  v.  (t)   given 
as 

9  s 
P 

and  A  is  the  logarithmic  amplifier  voltage  constant  determined 
approximately  by 
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»     +5000+R      (26xl„-3 

Aas( ) -^ ^  Volts  (3.13) 

where      R~      is    a   temperature    dependant    resistor  within   the    logarithmic 
amp  1 if  ier  modu le . 

I,       is    the    log   amplifier    reference    current    determined 
primarily   by      R    . 

The    concept   by  which   the    circuit   works    is    as    follows: 

Using   the    ion   chamber    current   as    the    input,    the    logarithm  of      i    (t) 

is    first    taken.      The    output   voltage      (v..(t))       is    a    linear   ramp 

f unct  ion  s  ince 

i  (t)  I       t/e 

vl(t)    =-Alog10^-=-Alog10(^e        P) 

k  k 

=  "A  logio  if  "  f  loHoe  (3'14) 

k    p 
v  (t)   is  then  differentiated  to  produce  a  d-c  voltage  V   where 


14 

Modular  Logarithmic  Amplifier  Model  4116    (Burr-Brown  Research 

Corp,  ,   1971)  ,   p. 3. 

A  more  complete  discription  of  the  logarithmic  amplifier  and 
its  parameters  is  given  in  Appendix  A. 

1  ft 

Since   the    input    current       (i    (t))       is    an   exponentially    increasing 

function  and   since    the    circuit    (Figure   3.4)    used   to  process    this 
function  has   a   limit    to   the    absolute  magnitude   of   the    input,    i    (t) 

can  only  remain   "on"   for  a   discrete   period  of   time      9    .      Therefore, 

in  order    to  be   absolutely   correct   when   defining      i    (t)      the   gate 

function  notation  should  be   used  where 


Ge  (t  "  9t/2}  =  l     °  *  l  * 


t 


Ge  (t  "  et/2}  =  °      t  <  0,  t  >  et 

is  the  time  of  the  input  pulse. 
Using  this  notation   i  (t)   becomes 
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v2  =  "R2ci  37  V0  =  "R2ci  3F("A  logio  t  •  f-  logioe> 

k         p 

R«C.A 
=  -^T1-  log„rte  (3.15) 

e  10 

p 

This   d-c  voltage    is    then   filtered   to  eliminate    any    low- 
frequency   noise   enhanced   by    the    differentiator.      The   resultant 

voltage     V  is    inversely   proportional   to   the   reactor   stable 

out 

period   and   approximately   proportional    to   the    reactivity   such   that 
R2C1A  R?C1A 

Vout  =  -T"  logioe  -  p  ~X~  logioe  (3-16) 

P  0 

As  previously  discussed  (Section  2.2)  the  relationship 
between   9   and   p   is  nonlinear.   Calibration  of  the  meter  may 
be  accomplished  by  producing  known  reactor  periods  and  marking 
the  meter  deflection  with  the  corresponding  values  of   p  via  the 
graph  in  Figure  3.6. 

A  comprehensive  analysis  of  the  circuit  design  is  given  in 
Appendices  A,  B,  C  and  D. 

3.4  Calibration 

The  output  voltage  in  the  circuit  of  Figure  3.4  is  a  d-c 
voltage  proportional  to  the  inverse  of   9   and  approximately 


c/e   r 


ls(t)  =  I0e   » 


Get(t-8e/2> 


If  the  gate  function  notation  is  carried  with   i  (t)   two  finite 
strength  impulses  will  appear  when  v1 (t)   is  differentiated  in  Equa- 
tion (3.15)  due  to  differentiation  of  the  gate  function. 

Although  these  finite  strength  impulses  are  of  importance  (partic- 
ularly since  they  show  up  on  the  output  when   i  (t)   is  switched  on 

and  off)  they  shall  not  be  carried  in  the  analysis „   It  should  be  under' 

stood,  however,  that   i  (t)   is  on  for  a  finite  period  of  time.   This 

s 

period  of   time      (9  )      is   discussed   further    in  Appendix  B. 
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(3.16) 


proportional  to   p   where 

R2C1A  R2C1A 

Vout  =  -~Y~   1OS10C  S  P  "I—  1O810G 
P  0 

As  previously  mentioned  the  relationship  between   Q   and 

P 

p      is   nonlinear.      Figure   3.6    illustrates    this    relationship   for   the 

Michigan   State   University  TRIGA  Mark  I    reactor. 

The   range   of   reactor   stable   periods    for  which    the  meter    is 

to  be    calibrated    is   10   sec   £9      £  40   sec.      This    range    corresponds 

P 

to  a  range  of  reactivity  between  40  and  19  cents  respectively. 
(F  igure  3  .  6) 

Calibration  of  the  reactivity  meter  is  accomplished  with 
the  apparatus  in  Figure  3.5. 


Reactor  Power 
Level  Indicator 
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Ion  Chamber 

Figure  3.5   Apparatus  used  to  calibrate  the  reactivity  meter 

The  reactor  is  first  made  critical  with  the  ion  chamber 
current  equal  to  I~.   A  control  rod  is  then  pulled  out  a  discrete 
distance  and  the  ion  chamber  current  becomes 

t/e 


1B(t)  =  V 


(1.1) 
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The  corresponding  d-c  voltage  on  the  reactivity  meter  will 
be 

RCA 

*«*--¥"  l°hoe  (3-16) 

P 
This  voltage  is  then  marked  on  the  meter. 

In  order  to  mark  the  reactivity  corresponding  to  V 

°      out 

the  graph  of  Figure  3.6  is  used.   As  this  graph  indicates  the 
reactor  stable  period  must  be  known.   This  value  is  determined  by 
using  the  reactor  power  level  indicator  and  a  timer  to  find  the 
power  level  doubling  time   9,.   The  period  can  be  calculated  from 
the  doubling  time  since 

p   t"  2 

This  process  is  repeated  until  the  range  of   p   is 
adequately  covered. 
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Figure   3.6     Graph   showing   the   relationship   between   reactivity   and 

reactor   stable   period    for    the  Michigan  State   University 
TRIGA  Mark  I   nuclear   reactor    (Courtesy   of  Gulf  General 
Atomic) 


CHAPTER  IV 
CONSTRUCTION  OF  THE  REACTIVITY  METER 

4.1  Circuit  Tests,  Performance  and  Specifications 

The  circuit  of  Figure  3.4  was  constructed  and  tested. 
Figure  4.1  shows  several  photographs  of  the  instrument. 

Testing  of  the  circuit  was  accomplished  by  connecting 
the  meter  input  to  the  ion  chamber  via  a  coaxial  cable.   The  re- 
actor was  then  made  critical  and  a  control  rod  withdrawn  to  pro- 
duce a  reactor  period  within  the  desired  range.   The  range  of 
the  meter  was  checked  for  consistancy  and  variation  due  to  noise. 
These  tests  revealed  the  meter  deflection  to  be  approximately  77 
percent  for  a  10  second  period  and  approximately  24  percent  for  a 
40  second  period  thus  adequately  covering  the  desired  range  of 

Q  .   The  maximum  variation  in  meter  reading  due  to  noise  was  less 
P 

than  +0.4  percent  of  full  scale. 

The  minimum  reactor  power  level  for  which  the  circuit  was 
designed  to  operate  is  ten  watts.   This  corresponds  to  an  ion 
chamber  current,  I„  =  3nA.   Values  of  input  current  below  this 
level  will  put  the  logarithmic  amplifier  in  saturation  and  there- 
fore cause  V     to  register  inaccurately  on  the  meter, 
out 

The  upper  power  level  was  chosen  to  be  200  watts  resulting 
in  an  ion  chamber  current  of  60  nA.  This  value  is  not  an  absolute 
upper  limit  and  the  meter  will  continue  to  read  accurately  up  to 
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a  reactor  power   level  of  1000  watts    (i    (t)    =  280  na)    at   which   time 

s 

the  negative  temperature  coefficient  inherent  in  the  reactor  will 
cause  the  meter  to  lose  accuracy. 

The  final  specifications  of  the  reactivity  meter  are  listed 
in  Table  4.1. 

Table  4.1   Specifications  of  the  reactivity  meter 


t/6 

Input:      Exponentially    increasing   current .». i    (t) 

■v  p 

.    10   seconds 

p  min 

Q                       

.    40   seconds 

p  max 

I                3   x 

-9 
10        Amperes 

s   min 

I                 60   x 

-9 
10        Amperes 

s  max 

(p)    in   cents 

rmm 

rmax 

Since  the  values  of   p   on  the  meter  are  calibrated  from 
a  graph,  the  only  error  will  be  that  variation  caused  by  noise  on 
the  output  voltage  of  the  filter  (Figure  3.4). 


CHAPTER  V 
CONCLUSIONS 

The  control  rod  calibration  curve  in  Figure  1.1  is  of  the 
type  for  which  the  reactivity  meter  was  designed  to  produced. 
Original  design  specifications  called  for  the  measurement  of  periods 
between  10  and  40  seconds  corresponding  to  reactivities  between 
40  and  19  cents.   As  the  graph  indicates,  these  values  cover  only 
a  small  portion  of  the  total  rod  calibration  curve. 

In  order  to  cover  the  entire  range  of  the  calibration 
curve  the  measured  values  of  rod  reactivity  must  be  added  in  incre- 
ments within  the  range  of  the  meter.   For  example:   Reference  to 
the  graph  of  Figure  1.1  shows  the  reactivity  corresponding  to  a 
rod  position  of  300  to  be  approximately  60  cents.   If  at  this 
position  the  reactor  were  made  critical  and  the  rod  pulled  out  to 
a  position  of  350,  a  reactivity  of  about  40  cents  would  be  recorded 
on  the  meter.   In  order  to  get  the  total  rod  worth  at  this  new 
position  this  value  of  reactivity  would  have  to  be  added  to  the 
previous  total  reactivity  of  60  cents  and  would  therefore  yield 
the  new  total  reactivity  of  approximately  1  dollar. 

It  should  also  be  noted  that  since  the  total  rod  worth  is 
determined  by  the  sum  of  incremental  reactivities  the  error  in 
measurement  of  each  incremental  rod  worth  becomes  very  important 
and  the  total  error  increases  with  each  sum. 

30 
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The  overall  results  of  the  tests  made  on  the  instrument 
reveal  that  the  reactivity  meter  works  well  for  the  range  of  input 
currents  and  reactor  periods  for  which  is  was  designed. 
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APPENDIX  A 
DESIGN  OF  THE  LOGARITHMIC    AMPLIFIER 


The    logarithmic  amplifier  used   in  constructing  the   circuit 
of  Figure   3.4   is   a  Burr-Brown  model  4116  modular   logarithmic 
amplifier.      Its   purpose    in   the   circuit    is    to  take   the    logarithm 
of   the    input    current      i   (t)      where 

t/e 
i.w  =  v     »  (i.D 

The   resultant   output    is   a   linear   ramp   function  given  by 

i    (t)  In      t/9 

vl(t)    =   -A  log,,,  -§—  -   -A  log1()     jtt  e        P 

"  "A  loPio  \  -  f  logioe  (3-14) 

k    p 
There  are  two  design  parameters  for  this  amplifier.   These 
are  the  voltage  scale  factor  A  and  the  reference  current   I  • 
Figure  A.l  shows  the  logarithmic  amplifier  with  exterior  connections. 
The  electrical  specifications  are  given  in  Table  A.l. 


The  design  equations,  graphs,  figures  and  specifications  in  this 
section  were  taken  from  the  pamphlet  Modular  Logarithmic  Amplifier 
Model  4116  (Burr -Brown  Research  Corp.,  1971),  pp.  1-5. 
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Table   A.l     Electrical   specifications    of   the   Burr-Brown  modular 

logarithmic  amplifier  model  4116    (current    input   only) 


Accuracy 

Accuracy,    percent   of   full   scale   with    current 

source    input    for   0.4nA  <.  I      <  400  (j,A +  1   percent 

Input 

Current   Source   Input    400   pA  to  400  ^,A 

Absolute  Maximum  Current    +10  mA 

Reference  Current   Range    (I,  ) 400  nA  £  I      <.  40  p,A 

K  K 

Voltage   Scale  Factor  Range    (A)     2/3   V   <.  A  <.  10  V 

Output 

Current    +5  mA 

Voltage    +10  V 

Output  Impedance   at   A  =  5    10  Ohms 

Stability 

Scale  Factor  Drift    (AA/°C)    +0.0005   A/°C 

Reference   Current   Drift    (&I    /°C) 

for  0.4  |j,A  <.  I,    <.  1  p,A   +  0.003   I   /°c 

Input  Offset  Current  Drift  (AI  /°C) . . .10  pA  at  25°C,  1  pA/°C 

Input  Noise  -  Current  Input  ....  1  pA  rms ,  10  Hz.  to  10k  Hz. 

Power  Supply 

Rated  Voltage  +  15  Vdc 

Supply  Drain  (Quiescent)  +  22  mA 
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■A/\ »+  15  V 


Current    Input       gl 


Voltage   Input         <> 
(optional) 


r 


3  + 

X     BB    -  4116 

1     Log   Amplifier 

2 

Y 


R 

— Vs 


T 


+  15  vi     d-c   c 


offset        6-  15  V 
adjust 


-0    V 


Figure   A.l     Burr-Brown  modular    logarithmic   amplifier  model 
4116  with   external   connections 


A   design    that    takes    advantage    of   the    full   range    of   out- 
put  voltage    is   best.      The   values    of     R        and     R        required    to 

x  y 

accomplish    this  may   be    obtained  by   first   solving  Equation    (3.14) 
for    the    full  +  10   volt    range   where 

I 


Vn  =    -A    log,„   -S^  =  10  Volts 

1  max  &10   I, 


(A.l) 


V,       .      =    -A    log1A  -^x   =   _10  Volts 
1  mm  10   I, 


(A. 2) 


result 


These    equations   may  be    added    to   solve    for      I        with    the 


I,     =    (I  I    .    ) 

k  max      rain 


(A. 3) 


The  corresponding  voltage  scale  factor  (A)  may  next  be 
found  by  substitution  of  Equation  (A. 3)  into  Equation  (A.l)  for 


I,  .   This  provides  an  expression  for   A  where 
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A  -  12 (A.4) 


In  Chapter  IV    it  was   established   that   the  value  of     I    . 
r  mm 

would   be      3nA      and      I  would  be      60nA.      Substitution  of    these 

max 

values    into   the   expressions    for      I,       and     A     yields 


I,    =  12.96  tiA  (A. 5) 

k 

A  =  15.75  Volts  (A.  6) 


At   this   point   a  problem    is    encountered.      Figure   A.  2    shows 
the   graphical   relationships   between      I,  ,   A      and    their   design   re- 
sistors.     Reference    to   these   graphs   and   to   the   specifications    in 
Table   A.l    indicate    that    the   values    of      I,       and     A      in  Equations 
(A. 5)    and    (A.  6)    are    outside    the   design   limits. 

A  compromise  must   be   made    to   reduce    the    output  voltage 

range    in   order   to  put      I,       and      A     within   their   respective    tolerances 

For   the    final   design   a  value    of     R      =  11  .3  MQ  was    selected.      This 

x 

corresponds    to   a  reference    current    of      .5  6jjlA. 

Using   this   value    for      I,       the   best  value    for      A      can   be 

calculated   by    letting    the   output   be   10  volts    for     I  .      Sub- 

J  °  r  s   min 

stitution    into  Equation    (A.l)    yields 

A  -  -10T  - =12 —  =  4.5  (A.7) 

i  min  -  3X10 

Reference    to   the    graph    in  Figure   A . 2b    indicates    a  value    of 
28   kfi   for     R        will  make      A     equal   to  4.5   volts. 

y 

The  actual  output  range  of  the  logarithmic  amplifier  may 
now  be  calculated  where 
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100  p,A 


10  p,A 


w*  1  M-A 


100  nA 


10  nA 


10  k    100  k     1  M     10  M     100  M 
R   (  ohms  ) 

y 


10  v 


8  V 


6  V 


4  V 


2  V 


(a) 


20  k 


40  k 


60  k 


80  k 


R   (ohms) 
x 


(b) 

Figure  A.  2   Design  curves  for  the  Burr -Brown  model  4116 

log  amplifier,  (a)   Reference  current  vs.   R  , 
(b)   Voltage  constant  vs .   R  y 
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V,  =   10  Volts  (A. 8) 

1  max 

I  -9 

,7  .    -  max  ,    c    -,  60x10 

=  4.36  Volts  (A. 9) 


The    design   parameters    of    the    logarithmic    amplifier  may 
be   summarized  where 

A      =  4.5   Volts 
R      =   28.    k  Ohms 

y 

I      =    .56  u.A 

R      =  11.3  M  Ohms 

x 


Input   Range    (i    (t))  3nA    to   60nA 

Output   Range    (v,(t))  10V    to   4.36V 
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APPENDIX  B 
DETERMINATION  OF  THE  RAMP  FREQUENCY  COMPONENTS 

The  ideal  output  of  the  logarithmic  amplifier  in  the  circuit 
of  Figure  3.4  is  a  linear  ramp  function  of  the  form 

vl(t)  =  -A  log10  ^  -  f  log10e  (3.14) 

k         p 

Differentiation   of   this    function  will  yield   a   d-c  voltage 

where 

v2  "  -R2cx  3T  Vc>  "  -R2ci  37  ("A  lo*io  if  -  f  logioe) 

k         p 
R2C  A 

-    e    lo«ioe  (3-15) 

P 

In  order  to  obtain  the  true  differential  as  in  Equation 
(3.15)  it  is  necessary  to  differentiate  all  frequencies  of  the 
ramp  function  v. (t) .   This  therefore  requires  the  knowledge  of 
the  range  of  ramp  frequency  components. 

The  highest  frequency  components  occur  when  the  ramp  func- 
tion is  on  for  the  shortest  period  of  time.   This  corresponds  to 
a  reactor  stable  period  of  ten  seconds. 

The  "on"  time  for  this  ramp  may  be  found  by  substituting 
the  limits  of   i  (t)   into  Equation  (1.1)  and  solving  for   9   where 


"t 

"max    min 


I    =  I  .  e  t      P  min  (1.1) 
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or 
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-9  -9      9t/l0 

60   X   10        =  3    X  10        e 


9      =  10£n   20  ^  30   seconds 


The  Fourier    series   half -range    cosine   expansion   for 

18 


f(x)    =  x      as    a   recurrent    triangular   wave    is' 


,.,   N  k       4k    .  rrx        1  3nx        1  5ttx 

f  (x)    =  x  -  j  "  "J  (cos  —  +  -  cos  —  +  —  cos  — 

TT 


+  ...+ l—  cos    (2n-l)rrX)        0<x<k 

(2n_1)  n-1,2,...  (B.l) 


Substituting   the    ramp    function   of  Equation    (3.14)    for 


f(x)      the   Fourier   series    expansion   for     v    (t)       is 


/    N  A     ■■  0      At    ,  .     n  0  t 

vl(t)    =   -A   log1()  £  —   log1()e   =    -A    log1()  -  -  - 


e-  A  logioe 


k        p  k  p  rain 


49tA  -  f  nt        1  3jTt      1  5nt 

— log1()e    cos   -  +  -  cos   —f  -  cos   — 

B  \  t  t  t 

p  mm  \ 

1  7rrt        1  9nt  1  lint 

+   -—    COS    — ! h    — -    COS    — 1 — "    COS 


49          et  8i          et     121            et 

l          1  (2n-l)TTt|            n        f        n          ln 

+  .  .  .+  2    cos 0  <   t   <   9      =30   sec. 

<2n-"  9t        /            a  -1,2,.!.                                        (B.2) 


The    radian    frequency    components    for    the  minimum  duration 
ramp   and   the   respective    amplitudes    of  each      (P.)      can  be   determined 


from  Equation    (B.2).      The  values    are 


18 

Samuel  M.    Selby,   CRC   Standard  Mathematical  Tables    (Cleveland, 

Ohio:      The  Chemical  Rubber  Co.,    1960),    p.    147. 
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46  A   log     e 


j.    =  —  =  0.1048   rad/sec,     P.    =  -~ —  =  7.9 

1       9-  1         2n 


t  TT    9 

p  min 

3rr  49^A    1Og10e 

uu~   =  —  =  0.314  rad/sec,      P»    -  — \ =^-  =  0.878 

9t  J       9n  9 

p  min 

c                                                          49  A  log     e 
uv    =  —  =  0.524  rad/sec,     Pc    =  — K —  =  0.316 

5        9t  5       25tt29        • 

p  mm 

_  49  A  log     e 

uo_    =  —  =  0.734   rad/sec,     Vn    =  — ~ «*-  =  0.161 

7        9t  7        49n29 


p  min 
1  lt>810' 


49  A  log     e 

uoQ  =  f11  =  0.943   rad/sec,     PQ  =  — ^ ^-  =  0.0975 

y        9t  9       81rr   9        - 

p  min 


Utt                                                       4e^A  lo8me 
ou   .    =  iJJL  =  1.15   rad/sec,     P        =  — ±— ^-  =  0.0653 

1     9t  U   121n  9    . 

p  min 

Since  the  amplitudes  of  higher  frequency  components  are 

less  than  1/144  of  the  amplitude  of  oj,   it  is  reasonable  to  omit 

them  from  the  analysis.   The  value  of  uo    is  therefore  the  highest 

radian  frequency  to  be  differentiated. 
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APPENDIX  C 
DESIGN   OF  THE   DIFFERENTIATOR 


Figure  C.l    is   a  magnitude   plot    for   the    real   and    imaginary 
components    of   an    ideal   differentiator  whose    transfer   function    is 


out 

V.    (10) 
in 


-  jKuo 


(C.l) 


V  ((A)) 

out 


UL) 

(b) 


Figure  C.l  Transfer  function  plot  for  an  ideal  differentiator, 
(a)   Imaginary  part,  (b)   Real  part 

A  differentiator  circuit  with  the  type  of  transfer  function 
in  Equation  (C.l)  has  the  disadvantage  of  enhancing  noise. 

If  the  range  of  frequencies  to  be  differentiated  is  limited, 
it  is  possible  to  reduce  unwanted  noise  by  using  a  compensated 
differentiator  whose  transfer  function  is 


19 

Albert   S.    Jackson,   Analog  Computation    (New  York:      McGraw-Hill 

Book  Co. ,    1960) ,   p.    147. 

42 


where 


and 


C.2 


Re 


Im 


V    (uu) 
out 

V.(uo) 
in 


out 
V.n(0)) 


Vout(^ 

v.n(.) 
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jKuj 


(jKu)  +  l)(jKuo  +  1) 


(C.2) 


o  2  2 
2K  uo 


2  2  2     2  2 
(1  -  KV)   +  4K  u) 


iKood  -  K2gj2) 

2  2  2     2  2 
(1  -  K  u)  )   +  4K  uo 


(C3) 


(C.4) 


These  real  and  imaginary  components  are  plotted  in  Figure 


Figure  C.2  Transfer  function  plot  for  a  compensated 

differentiator,  (a)   Imaginary  part,  (b)   Real  part 


In  Appendix  B  it  was  determined  that  the  maximum  radian 
frequency  to  be  differentiated  in  the  linear  ramp  function  was 
uo.  1  =  1.15  rad/sec.   It  is  therefore  desirable  to  make   1 /K  greater 


than  1 .15  rad/sec, 


Figure  C.3  is  a  circuit  diagram  of  the  differentiator 


20 


20 


RCA  Linear   Integrated  Circuits    (RCA,    1967),    pp.    75-78 
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*-AA-|( 


Figure   C.3     Equivalent    circuit    of   a   compensated   differentiator 


The    transfer    function  may  be    derived   as    follows 
Define  : 


R 


Zr  =  R„//C„    = 


f  2        2        juuR9C     +  1 


2    2 


Z      =  R.    + 


r  1        juoC 


(C5) 
(C6) 


V    ,   V.      and     V    „      may   be   written  as 
e'      i  01 


(Z      +  z      )//Z.  V      (Z    //Z.) 

v     _  v  f         0l  i  ,      01     r        r 

e   "      in  Z      +    (Zc  +  Z„.)//Z.  Z,  +  Z^  .   +  Z    //Z  . 

r             f           Oi           i  f           Oi           ri 


(C7) 


V.    =  0 

l 


V01   =    -V*)(Ve    -V 


(C.8) 
(C.9) 


Substitution   of     V        and     V.      into   the   expression   for 

e  l 


'oi     yields 


V 


-V.nA0(.)(Zf+ZQ.)Z, 


01         (Zf  +  Zo.)Z.   +  Zr(Zf  +  Z0.    4-  Z.)    +  A0(oJ)Z.Zr 


(CIO) 


If      R     -»  oo      the   exact   expression   for    the    transfer   function 


can  be   written  as 
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V  Zrt.Z.  -  A.(u))Z  .Z. 

out  _  Oi  1    0V  '    1  f (      ... 

V.    "  A-CuOZ.Z  +  (Zc   +  Z^.)Z.  +  Z  (Z.  +  Z  .  +  Z.)         ^    ; 
in     Uv  '  i  r    v  f    Oi   i    rv  f    Oi    i 

This  expression  may  be  simplified  if  the  following  assump- 
t  ions  are  made : 
Assume : 

Z  .  »  Z  //Z  (C12) 

i     f   r 

ZQ.  «  Zf  (C.13) 

Zf  «  A0(uo)  (C.14) 

Zr  «  AQ(UJ)  (C.15) 


21 
Equation  (C .  9)  may  then  be  written 


V  -z 

out  _  f 

V.  Z 

in  r 


(C.16) 


When  the  expressions  for  Z    and  Z..   are  substituted  into 

r  r        f 

this  equation  the  differentiator  transfer  function  becomes 


-21it   =  2   1 

Vin  O^i00  +  DU^S*0  +  1} 

If      R-jC..    =  R<>C„      then   the   magnitude   plot    in  Figure  C.2 
will   correspond   to  Equation    (C.17)    with     1/K  =  1/R.C    . 

For  a  differentiator  error   of  one   percent  maximum   it    is 
necessary   to  establish  a  corner  frequency  about   ten  times   the 
maximum  frequency  to  be   differentiated.      In   other  words 

=  10  to  =  10  w  (C.18) 


R1C1       R2C2 


21 

Ibid. ,   pp.    75-78. 
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22 
where      uo,  ■,    =  1.15    rad/sec      and  was    determined    in  Appendix  B. 

A  trade    off   of   accuracy   vs.    noise    is   encountered    if  Equa- 
tion   (C.18)    is   used   to   design   the   differentiator.      The   prupose    of 
the    low-pass    filter    in  Figure    3.4    is    to  attenuate    the    noise   passed 
by   the    differentiator.      It   was    discovered    that    for     1/R-C      =  10  m 
=  11.5    rad/sec      the    filter   could   not   effectively   reduce    the   noise 
without    causing   aneccessive   delay    in   the   d-c   response   due    to    its 
long   time    constant.      A  decision  was   made    to    limit    the    filter    time 
constant   to  0.6  second.      This  meant   the   corner   frequency   of  the 
differentiator  had   to  be   reduced   until    the   noise    level  was    tolerable. 
The   experimental  resultant   corner  frequency  was      1/R-C      =3.8   rad/sec 
or  about   three   times   as   great   as    the  highest    frequency   to  be   dif- 
ferentiated. 

The  magnitude    of   the    error   caused  by    inexact    differentiation 
of   the  maximum  ramp   frequency  will  be    small   since    the   amplitude    of 
this    frequency   component    is   only  1/121   of   the  major  ramp   frequency. 


22 

Jackson,    Analog  Computation,    p.    147. 
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APPENDIX  D 
DESIGN   OF   THE    FILTER 

Figure   D.l    is   a   schematic   of   the    output    filter    in  Figure 
3.4. 

The    filter    is    a   simple   RC    low-pass    filter  with   a  unity 
gain   output   buffer   stage    to  prevent    loading. 


out 


d-c   offset 

adjust 
15  V 


Figure   D.l     Output    filter   for   the    reactivity  meter 


The    transfer   function   for    the    filter    is 


V  fa) 

out 


V.n(u>)  jRfCfu)  +   1 


(D.l) 


Since   the   output   voltage    of   the    differentiator  (V„)       is 

a   d-c  voltage   with   superimposed    low-frequency  noise    it  is   desirable 

to  make   the    filter    time    constant    long  enough   to   reduce  this    noise 
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while  at  the  same  time  have  a  reasonable  response  time  to  the  d-c 
voltage . 

Final  tests  of  the  reactivity  meter  indicated  that  a  filter 
time  constant  of  RfC   =0.6  seconds  reduced  the  output  noise  level 
on  the  meter  to  less  than  +0.4  percent  full  scale.   This  value 
for  RfCf   also  gave  a  d-c  response  time  of  about  3  seconds. 
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APPENDIX  E 
SPECIFICATIONS   OF   THE   COMPENSATED   ION   CHAMBER 

The  TRIGA  Mark  I   nuclear    reactor   at   Michigan  State   University 

contains    a  Westinghouse   type   WL-8501    compensated    ion   chamber.      Table 

23 
E.l   tabulates    the   specifications    of   this    chamber. 


Table  E.l      Specifications    of    the  Westinghouse    type   WL-8501 
Compensated    ion  chamber 


Maximum  Voltage   Between  E  lectrodes    (dc) 1000  V 

Operating  Voltage    for  Boron  Chamber    (dc) 580  V 

Operating  Voltage    for  Gamma  Chamber    (dc) 0   to   -37  V 

13  n. 
Output   Resistance,    Signal   to  Case,  Minimum 10        Ohms 

Output   Capacitance,    Signal   to  Case,    Approximately 130   pF 

2  10  2 

Thermal  Neutron  Flux  Range 2.5x10      to  1.5x10        n/cm   /sec 

Thermal  Neutron  Sensitivity 1  .5x10  A/nV 

-13 
Gamma  Sensitivity,   Compensated 2.5x10  A/R/hr 


Since    the    output    resistance    of    the    ion   chamber    is    greater 

13 
than  10        Ohms    the    chamber  may   be    treated   as    a   current    source. 


23 

Westinghouse   WL-8501    (New  York:      Westinghouse  Electric  Co., 

1965)  ,    pp.    1-2. 
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